Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 5 February 2008 (MN WF$t style file v2.2) 



Detached shells as tracers of AGB-ISM bow shocks 



C. J. Wareing 1 *, Albert A. Zijlstra 1 , Angela K. Speck 2 , T. J. O'Brien 1 , 
Toshiya Ueta 3 , M. Elitzur 4 , R. D. Gehrz 5 , F. Herwig 6 , H. Izumiura 7 , 
. . . M. Matsuura 8 , M. Meixner 9 , R. E. Stencel 10 , R. Szczerba 11 

1 Jodrell Bank Centre for Astrophysics, University of Manchester, Oxford Road, Manchester, M13 9PL, UK 
2 Department of Physics and Astronomy, 223 Physics Building, University of Missouri, Columbia, MO 65211, USA 
O ! 3 NASA Ames Research Center/USRA SOFIA Office, Mail Stop 211-3, Moffett Field, CA 94035, USA 
O^l ■ "^Physics & Astronomy Department, University of Kentucky, Lexington, KY 40506, USA 

t— H ' Department of Astronomy, School of Physics and Astronomy, 116 Church Street, S.E., University of Minnesota, 
Minneapolis, MN 55455, USA 

6 Theoretical Astrophysics Group, LANL, Los Alamos, NM 81545, USA 

7 Okayama Astrophysical Observatory, National Astronomical Observatory of Japan, Kamogata, Asakuchi, 
Okayama 719-0232, Japan 

^Division of Optical and IR Astronomy, National Astronomical Observatory of Japan, Osawa 2-21-1, Mitaka, 

Tokyo 181-8588, Japan 
9 STScI, 3700 San Martin Dr., Baltimore, MD 21218, USA 

10 Department of Physics and Astronomy, University of Denver, Denver, CO 80208, USA 
11 N. Copernicus Astronomical Centre, Rabianska 8, 87-100, Torun, Poland 



submitted to MNRAS Letters 



ABSTRACT 

New Spitzer imaging observations have revealed the structure around the Mira 
variable star R Hya to be a one-sided parabolic arc 100 arcsec to the West stretching 
from North to South. We successfully model R Hya and its surroundings in terms of an 
interaction of the stellar wind from an asymptotic giant branch (AGB) star with the 
interstellar medium (ISM) the star moves through. Our three-dimensional hydrody- 
namic simulation reproduces the structure as a bow shock into the oncoming ISM. We 
propose this as another explanation of detached shells around such stars which should 
be considered alongside current theories of internal origin. The simulation predicts 
the existence of a tail of ram-pressure-stripped AGB material stretching downstream. 
Indications for such a tail behind R Hya are seen in IRAS maps. 

Key words: stars: AGB and post-AGB - stars: individual: R Hya - stars: mass-loss 
- ISM: structure. 



1 INTRODUCTION 

Large detached shells have been observed around several 
asymptotic giant branch (AGB) star s. They have been 
seen in IRAS images of dust emission dWaters et alJll994t 
[izvirrm jja^^t,^ |l997|), CO line emission lojofesonet^d] 
Il996l: ISchoier. Lindavist fc Olofssonl 120051) . and in a few 
cases [Na i] and [K i] emission as well as in th e optical contin- 
uum jGonzalez Delgado et aLllioOll 120031) . lOlofsson etafl 
J1990D suggested that such shells are the result of mass-loss 
variations and in particular, a thermal pulse or He-flash. 
During a He-flash, an intense short-lived mass ejection is 
driven by the star reaching a critical luminosity. Thermal 
pulses are separated by phases of quiescent hydrogen burn- 
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ing lasti ng 10 4 -10 5 yr. The s t ellar evolution tracks calcu- 
lated by IVassiliadis fc Wood! (Il993ft confirmed that mass- 
loss fluctuations during the thermal pulse cycle can lead to 
detached circumstellar shells. Hydrodynamic simulations by 
ISteffen fc Schonbernerl i2000t) showed that a brief period of 
high mass loss can translate into a geometrically thin shell 
expanding around the star. This has become the standard 
explanation of detached dust shells aro und stars and obser- 
vations have been interpreted as such JZiilstra et alJll992l : 
ISpeck. Meixner fc Knappll2000lL 

A separate explanation for large detached shells 
is the interaction of t he AGB wind with the inter - 
stellar medium (ISM ) lYoung. Phillips fc Knapd Il993l) . 
IZiilstra fc Weinberger! (120021) used this for a giant (~ 4 
pc at D ~ 700 pc) detached shell surrounding an M3 
III AGB star. They proposed that its AGB wind has 
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been stopped by the surrounding ISM and the swept-up 
'wall' is now expanding at the local sound speed. Simu- 
lations by Villavcr. Garcia- Segura fc Manchadol (|2003h and 
IWareing et alJ i2006f) confirmed the viability of this mecha- 
nism. It is likely that both mechanisms occur, but whether 
the external mechanism of an ISM wall, or the internal 
mechanism of long-term mass-loss variations is the domi- 
nant cause of detached shells is not known. Schoier et al. 
(2005) found that the derived masses of the shells increase 
and the expansion velocities decrease with increasing radial 
distance from the star. They suggest that the shell is sweep- 
ing up surrounding material from an earlier mass loss phase. 
However, ISM sweep-up could yield similar effects. 

A difference between the two mechanisms is that the in- 
ternal one will normally give a spherical shell, while external 
mechanism will give a shape which depends on the motion 
of the star through the ISM. This is a testable prediction, if 
the proper motion of the star is known. Here we show that 
the detached shell around the Mira variable R Hya is due to 
an ISM bow shock. 



2 OBSERVATIONS 

R Hya is one of the b rightest Mira variables on the sky. 
Hash imoto et al.1 il99ST> found evidence in IRAS data for 
a detached shell 1-2 arcmin from the star. The star itself 
coincides with an IRAS point source due to an inner dust 
shell. This inner shell is also detached from the star, with 
an inner radius of about 60 stellar radii. Since its discovery 
in AD 1662, the pulsation period has decreased from 495 
days to 385 days, a ttributed to non-linear pulsation or a 
recent thermal pulse dZiilstra. Bedding fc Mattel2002T) . The 
inner detachment suggests a mass-loss interruption occurred 
around AD 1800. The outer shell has a dynamical age of 
around 8000 yr. This small time difference causes problems 
if one assumes that both detached shells are due to thermal 
pulses (Zijlstra et al. 2002). 

The reported proper motion from the Hipparcos cat- 
alogue is 60.73 masyr " 1 West and 11.01 masyr" 1 North 
PCTrjrnan_et alj ll997l) . The radial velocity is —10.4 kms -1 
( Wilson! I1953T) . The distance to R Hya is un certain, asdis- 
cussed in section 4 of Zijlstra et al. (2002). lEggenl (1985) 
argues for a distance of 165 pc based on a proper motion 
companion, but the data this is based on appears to be un- 
published. The Hipparcos non-detection favours a larger dis- 
tance. We have adopted a distance of 165 pc. This distance 
puts R Hy a just beyond the edge of the Local Bubble in this 
direction llLallement et alJl2003f) . At our chosen distance, 
the proper motion is equivalent to a transverse velocity of 
48.5 kms -1 , and a space velocity of 49.5 kms - . This is 
within 2a of the average space v elocity for AGB stars, 30 
kms -1 jFeast fc Whitelockll200ol) . 

High angular resolution observations were ta ken with 
the Spitzer Space Teles cope iWern er et all I2004T) as part 
of the MIRIAD project JUeta et alJl200eT) . The Multiband 
Imaging Photometer aboard Spitzer (MIPS; Rieke et al. 
2004) image is shown in Fig. The central unobserved re- 
gion in the images avoids the central star which would satu- 
rate the detec tors. For details of the data reduction process, 
please consult lUeta et alJ J2006T) . 

These images show that the material around the cen- 



tral star is n ot a circular shell as su ggested by the IRAS 
observations llHashimoto et aDll998l) . but instead consists 
of a parabolic arc to the West of the star and regions of 
emission to the East. The arc appears brightest in the direc- 
tion of the proper motion, fading as it stretches away North 
and South. At its brightest point, the arc is 100 arcsec from 
the central star. The arc is 300 arcsec wide across the star. 
There appears to be no emission further West of the arc, 
suggesting there is no material related to R Hya West of 
the arc. It is not clear from these images that the arc is 
detached from the star as the region between the arc and 
star is in the central unobserved region. To the East of the 
star, several low-brightness regions of emission are present 
up to 300 arcsec from the star. The location and brightness 
of these regions would suggest they are associated with the 
star. Our simulations closely reproduce the morphology of 
the arc and the Eastern emission. 



3 SIMULATIONS 

Our simulations have been performed using a parallelised 
computational fluid dynamics program designed to solve the 
standard Euler equations of hy drodynamics using a second- 
order Godunov scheme due to iFalld (|l99ljj) . It is in three 
dimensions using Cartesian coordinates and includes the ef- 
fect of radiative cooli ng above 10 4 K via a paramet ric fit 
to the cooling curve of lRavmond. Cox fc Smith! jl976). The 
timestep is defined by the Courant-Friedrichs-Lewy condi- 
tion with a Courant number of 0.5. We also include a nu- 
merical viscosity to avoid on-axis issues in the simulations. 

Using this numerical scheme, we have modelled the in- 
teraction of a stellar wind ejected from a star as it moves 
through the ISM. The simulation is performed in the frame 
of reference of the star and the motion through the ISM is 
like that of an oncoming wind - hence our 'tw o-wind' model- 
Furth er details of the model are explained in I Wareing et all 
(l200rjh . We model a period of 5.5 x 10 4 years on the AGB 
and have used a simulation grid of 200 3 cubic cells (cell size 
of 4 x 10~ 3 pc) producing a grid of 0.8 pc in each direction. 

We set the AGB wind parameters of the model with 
a mass-loss rate of 3 x 10 -7 Mq yr -1 and a velocity of 10 
kms -1 (Zijlstra et al. 2002). The temperature of the wind is 
set at 10 4 K, which i s the lowest va lue for which the cooling 
function is defined (IWareind 12001 . The real temperature 
will be considerably less than this and indeed the tempera- 
ture of the undisturbed AGB wind in the simulation is on 
the order of a few tens of K. Our simulations do not model 
dust physics or radiation transport. The gas pressure in both 
winds is calculated assuming an ideal gas equation of state 
with an adiabatic index of 5/3 in both winds, effectively 
ignoring molecules. 

The position of the bow shock 0.08 pc ahead of the star 
at our adopted distance can be understood in terms of a ram 
pressure balance. This balance predicts a local ISM density 
of 0.6 Hem -3 which we take as the ISM density in our sim- 
ulations. The ISM is modelled as a warm neutral medium 
with a temperature of 8 x 10 3 K. The stellar motion corre- 
sponds to a Mach number of 3.65 and the AGB wind has a 
Mach number of 0.74. At 165 pc, R Hya is located 101 pc 
above the Galactic plane. ISM density, a nd hence pressure, is 
constant at jih = 2 cm" " 3 JSpitzerll978h up to a scale height 
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Figure 1. The background-subtracted, mosaicked MIPS colour maps of R Hya at 70 /im (left) and 160 /an (right) JUeta et alj f2006V 
The images are zeroed at the position of the star (indicated by the 'star') with North up and East to the left. The proper motion is 
indicated by the arrow in the bottom left. The parabolic arc is marked for clarity by the dashed line. Colour-scaling of surface brightness 
(in Mjysr - 1 ) is provided by the colour bar. 



of 100 pc feinnev fc Merrifieldll998l) . Above this, density 
drops off exponentially, giving an ISM density in the region 
of R Hya of 0.74 Hcm~ 3 in accidental agreement with our 
simulation. This is als o consistent with t he estimate of 0.4 
Hem -3 calculated bv lUeta et alJ i2006l) . Densities within 
the Local Bubble would be very much lower: the bow shock 
supports a location of R Hya beyond this region. 



4 RESULTS 

Fig.|!2]shows density on the plane in which the star is located 
40 000 years into the AGB phase. This point in time is an 
arbitrary choice since the wind structure reached a stable 
state after 25 000 years into the simulation. The image is at 
9.5° to the actual line of sight to R Hya. 

The wind from the star has driven a shock into the ISM. 
In the case of a stationary star, this is a spherical shell of 
AGB wind material sweeping up shocked ISM material. In 
the moving case, the shock forms into a bow shock expanding 
ahead of the star. Eventually, the bow shock reaches a max- 
imum distance ahead of the star (after 25 000 years in our 
simulation) which can be understood in terms of a ram pres- 
sure balance. Material is being ram-pressure-stripped from 
the head of the bow shock into a tail stretching downstream. 
The 'bow shock and tail' structure is the same as we have 
seen in other simulations considering a range of space veloc- 
ities and wind parameters. We note this as support for the 
formation of this general structure as a convergent result. 

The width of the bow shock across the star is approx- 
imately 0.24 pc agreeing with the observations. If the bow 



shock is a strong shock, the temperature of the shocked ma- 
terial at the head can be predicted by T ~ 3/16 mv 2 /k — 
33 000 K, where m is taken as the average gas component 
mass of TO x 10 -27 kg, v is the speed of the central star 
relative to the ISM and k is the Boltzmann constant. The 
temperature in the simulation at the head of the bow shock 
is found to be in agreement with this prediction. The sur- 
vival of dust above temperatures of approximately 1500 K 
is questionable, but at the low densities in the simulation, 
the dust and gas are decoupled and evidently the dust tem- 
perature can be much lower. 



5 DISCUSSION 
5.1 The bow shock 

The two-wind model can reproduce the appearance of the 
circumstellar structure around R Hya, with physical dimen- 
sions matching that of the bow shock. The head of the bow 
shock is in good agreement with the direction of motion. 
Small regions of emission downstream can be explained in 
terms of regions of higher density ISM encountered by the 
stellar wind. Such high density regions in the I SM have been 
show n to survive ablation by stellar winds iPittard et alJ 
120051) . 

In our simulation, we find a high temperature of 35 000 
K at the head of the bow shock, with material cooling 
rapidly as it moves down the tail. The high te mperatures are 
consi s tent with the ob served Ha emission iGaustad et alJ 
120011 lUeta et alJl20o3) . The high temperatures will affect 
the dust by collisional heating. However, it is impossible for 
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Figure 2. An image showing the gas density in the plane in which 
the star is located, parallel to the direction of motion, 40 000 years 
into the simulation. The ISM is flowing in from the right. This 
image is at 9.5° to the actual line of sight to R Hya. The colour 
bar is in units of 10 -6 Mq pc — 3 where 2.4 X 10 5 corresponds to 
10 Hem -3 . The image is 0.8 pc on a side, which corresponds to 
the physical image scale in Fig. Q at 165 pc. 



us to say whether this heat from the local gas is more impor- 
tant than the stellar radiation or the interstellar radiation 
field (Speck et al. 2000) in controlling the dust temperature. 

The Spitzer and IRAS images indicate that the bow 
shock region is detected through dust emission; the fact that 
the emission is strongest at 60/im indicates a dust tempera- 
ture of order 60 K. But other emission processes may play a 
role. Shocked regions can sho w [O il 63/xm and 146/xm lines. 
The IRAS lOO^rn detection Jrlashimoto et al.lll993) . in a 
band without strong lines, shows that dust contributes to 
the far-IR emission. 



5.2 Other cases: mass-loss variations versus ISM 
interaction 

Cases of detached shells around AGB stars include U Ant 
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U Hya jWatersetal.lll994h. YCVn 
t|), R Hya jHashimoto et al.lll99SI) and 
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IRAS02091+6333 ijZiilstra fc Weinbergerll2002f ) observed in 
the far-infrared; detached molecular (CO) gas shells around 
TT Cyg, S Set, R Scl, U Ant, U Cam, V644 Sco and DR 
Ser (Olofsson et al. 1996, Schoier et al. 2005). All are carbon 
stars, apart from R Hya and IRAS02091+6333. It has been 
suggested that the thermal pulse scenario may only lead 
to mass-loss spikes for carbon stars (Schoier et al. 2005) 
but this is very much dependent on the assumed mass loss 
prescription which is not well known. IRAS colours indicate 
that detached shells do also exist around oxygen-rich stars 
( Ziilst ra et aljll992T) and many more stars are likely to show 
large shells (Young et al. 1993). 

An interesting comparison can be made between R Hya 
and U Hya. U Hya has a thin cir cular shell with a radius 
of 120 arcsec iWaters et al]ll994h . It has a space velocity 



of 50 kms~ and the parallax gives a distance of 162 pc 
iPerrvman et alJll997l) . The circular nature of the shell sug- 
gests an internal origin such as the m a ss-los s variations pro- 
posed by, for example, IZiilstra et al] 1119921) and Schoier et 
al. (2005) . The age of the shell in this particular case is 6000 
years. The ISM interaction may be radially further away 
from the star. 

The star TT Cyg is surrounded by a ci rcular detached 
shell at a radius of approximately 40 arcsec llOlofsson et al] 
1996); the physical size of the shell is very similar to U Hya, 
in view of the larger distance of TT Cyg. The thin symmet- 
rical appearance supports an internal origin. Interestingly, 
the space motion of TT Cyg (50 kms -1 ) is almost all in the 
radial direction away from us and any bow shock formed by 
an interaction with the ISM would appear circular on the 
sky. But the fact that the slight offset of the star from the 
centre of the shell is at right angles to the direction of proper 
motion, favours an internal origin in mass-loss variations. 

For other stars there is insufficient data to decide on the 
cause of the detached shells. It may be that the well-studied 
carbon stars with detached shells are mostly due to inter- 
nal mechanisms, i.e. thermal pulses, while the fainter, less 
studied shells are dominated by ISM interactions. CO obser- 
vations of shells around carbon stars reveal relatively high 
velocities (about 20 kms -1 ) which favours thermal pulse 
origins (Olofsson et al. 1996, Schoier et al. 2005). Multiple 
detached shells d o exist in the AGB phase of ev olution, e.g. 
S Set and U Ant ( Gonzal ez Deleado et al 1 l2003l) and also R 
Hya: thus, both mechanisms may occur simultaneously. 

Young et al. (1993) in their analysis of 76 AGB stars 
resolved by IRAS found no evidence for distortions by inter- 
action with the ISM. Since R Hya is included in their sample, 
this lack of evidence can be attributed to the poor spatial 
resolution and image quality of the IRAS instrument. 

R Hya can be considered a typical case for mass-loss 
rate and ISM density although its space velocity is perhaps 
high. Lower velocity objects will have a bow shock and with 
higher mass-loss rates this will be located further from the 
star. In a zero velocity case, the bow shock transforms into a 
spherical swept-up shell of ISM material (Young et al. 1993, 
Speck et al. 2000, Zijlstra & Weinberger 2002). We predict 
that all AGB stars will show some degree of an AGB-ISM 
interaction, although bow shocks will be rarer. The interpre- 
tation of a detached shell in terms of a mass-loss variation 
must be considered with this ISM interaction in mind. 

5.3 Mass loss history 

This model has shown that information usually gleaned from 
circumstellar dust shells around Mira variables can no longer 
be inferred in this situation. The bow shock has destroyed 
any mass loss history older than about 10 4 years in the case 
of R Hya. Higher mass- loss rate and/or lower lower ISM 
density could increase this timescale as found by Young et 
al. (1993) and Zijlstra & Weinberger (2002). 

The simulations predict the occurence of a tail, consist- 
ing of swept-back ISM and stellar wind gas. The mass loss 
history can in principle still be traced down the length of 
the tail. IRAS maps provide some indication for a tenuous 
detection of material downstream of R Hya as shown in Fig. 
[3] This material stretches up to 30 arcmin away or 1.5 pc at 
a distance of 165 pc. At 49.5 kms -1 this implies a minimum 
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Figure 3. A figure showing an IRAS 60fim observation of the 
area around R Hya. North is up and East to the left. Evidence 
for a tail of material is indicated by the ellipse. 



tail age of 30 000 years. Adding the 25 000 years it takes to 
form the stable bow shock, to represent the travel time from 
the star to the bow shock and down the tails, we predict R 
Hya has been losing mass for at least 55 000 years. 

If we consider it has taken 25 000 years to form the bow 
shock and after this the bow shock is in a steady state, an 
appreciable amount of mass is in the bow shock. The esti- 
mate from the simulation is 2.6 x 10~ 3 Mq in the region 
of the bow shock defined as a hemispherical upwind shell 
centred on the star with inner radius of 0.08 pc and a thick- 
ness of 0.1 pc. Our analytical estimate of the mass in the 
bow shock is consistent with this estimate, suggesting there 
is four times more stellar mass in the bow shock than ISM 
mass. At a dust temperature of 60 K and assuming a dust to 
gas ratio of 100, we predict a 100 /im fl ux of 15 Jy, in bro ad 
agreement with the fluxes observed bv lUeta et al J J2006T) . 

We use a constant stellar mass-loss rate and wind veloc- 
ity. However, the dust shell models of various people (Zijl- 
stra & Weinberger 1992, Schoier et al. 2005) indicate that 
these quantities are variable, particularly during the ther- 
mal pulses. Models indicate that in the course of thermal 
pulses the stellar radius changes temporarily by a factor of 
2, and the mass-loss rate by an order of magnitude, causing 
a density and velocity spike moving radially away from the 
star - considered as the origin of detached dust shells. Such 
events could severely disturb the pressure balance that dic- 
tates the well defined location of the bow shock when one of 
these 'wind'-shells runs into the bow shock. There is likely 
to be a complex time dependence of the wind-ISM interac- 
tion on the timescales of the thermal pulses, as has been 
considered in 2D simulations by Villaver et al. (2003). 
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